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Abstract  —  This  work  provides  a  design  method  to 
achieve  a  specified  probability  of  resolution  for  two  tar¬ 
get  localization  via  a  wireless  network  of  proximity  sen¬ 
sors.  Proximity  sensors  simply  report  a  single  binary 
value  indicating  whether  or  not  a  target  is  near.  The 
design  provides  the  density  and  threshold  settings  to 
achieve  the  given  probability  of  resolution  for  targets 
separated  by  a  specified  distance.  Simulations  are  in¬ 
cluded  that  demonstrate  that  at  the  designed  sensor  den¬ 
sity  and  threshold  values,  the  actual  percentage  of  tar¬ 
gets  resolved  achieves  the  desirable  level  of  resolution 
for  moderate  to  large  target  separations. 

Keywords:  wireless  sensor  networks,  source  localiza¬ 
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1  Introduction 

Proximity  (or  binary)  sensors  are  simple  devices  that 
report  whether  or  not  they  detect  an  object.  A  binary 
sensor  can  consist  of  a  passive  device  such  as  an  acous¬ 
tic  microphone  where  it  reports  a  detection  when  the 
energy  of  the  signal  computed  over  an  integration  pe¬ 
riod  exceeds  a  threshold.  In  essence,  it  acts  as  a  tripwire 
sensor.  These  sensors  have  proved  useful  in  a  number  of 
applications  where  they  cue  more  expensive  devices  to 
interrogate  the  scene  using  a  more  expensive  and  infor¬ 
mative  sensor  upon  detection.  This  paper  investigates 
the  feasibility  of  a  network  of  proximity  sensors  to  local¬ 
ize  targets.  In  such  a  network,  the  no  detection  report 
is  as  valuable  as  a  detection  report.  Previous  work  has 
revealed  the  potential  of  target  localization  via  a  mesh 
network  of  proximity  sensors  by  proposing  localization 
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methods  for  the  single  target  case  [1,  2,  3,  4,  5,  6].  Re¬ 
cent  work  presents  maximum  likelihood  (ML)  localiza¬ 
tion  algorithms  for  multi-target  localization  when  the 
number  of  targets  is  known  a  priori  [7].  When  the 
number  of  targets  is  not  known  a  priori,  the  finite-set 
based  Probability  Hypothesis  Density  (PHD)  filter  [8], 
which  can  provide  reasonable  estimates  of  the  number 
of  targets  and  the  target  states,  might  prove  useful  for 
proximity  sensors.  To  understand  the  feasibility  of  any 
proximity-based  localization  method,  one  needs  to  de¬ 
velop  fundamental  performance  bounds  that  determine 
how  well  a  network  of  proximity  sensors  can  localize 
targets  given  the  operation  parameters  of  the  sensor 
network.  Here,  the  operation  parameters  include  the 
density  of  the  sensor  nodes  and  the  threshold  settings 
for  the  sensors. 

To  localize  targets  effectively,  one  must  be  able  to 
resolve  them.  To  the  best  of  our  knowledge,  the  no¬ 
tion  of  target  resolution  has  not  been  investigated  for 
a  proximity  sensor  network.  This  paper  presents  our 
initial  work  to  understand  target  resolution  for  such 
networks,  and  to  determine  the  operation  parameters 
to  achieve  resolution  of  targets  separated  by  a  specified 
distance.  To  this  end,  the  work  exploits  a  physically 
justifiable  sensor  model  that  determines  the  probabil¬ 
ity  of  detecting  targets  based  upon  the  location  of  the 
targets  relative  to  each  sensor.  As  an  initial  investi¬ 
gation,  a  number  of  simplifying  assumptions  are  made 
such  as  the  location  of  all  sensor  nodes  are  known  ex¬ 
actly,  the  power  of  radiation  from  all  target  sources 
are  equal  and  known,  the  propagation  loss  exponent 
is  known,  no  communications  errors,  perfect  time  syn¬ 
chronization,  etc.  Once  the  fundamental  performance 
bounds  for  the  simplified  case  is  developed,  the  assump¬ 
tions  can  be  relaxed.  Finally,  the  analysis  in  the  paper 
assumes  that  the  mesh  network  is  developed  by  tossing 
out  the  sensor  nodes  randomly  so  that  they  distribute 
via  a  statistical  uniform  distribution  over  the  coverage 
area.  Other  work  can  investigate  performance  bounds 
when  the  sensors  are  intentionally  placed  to  achieve  op- 
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timal  performance. 

The  paper  is  as  organized  as  follows.  Section  2  pro¬ 
vides  details  about  the  sensor  model,  and  Section  3  dis¬ 
cusses  the  ML  localization  method  and  demonstrates 
how  the  likelihood  surfaces  can  resolve  two  targets  un¬ 
der  various  conditions.  Section  4  uses  the  insight  from 
the  previous  section  to  devise  a  strategy  to  design  the 
operational  parameters  so  that  the  desired  level  of  tar¬ 
get  resolution  can  be  achieved.  The  effectiveness  of  the 
design  strategy  is  investigated  via  simulations  in  Sec¬ 
tion  5.  Finally,  Section  6  provides  concluding  remarks. 

2  Sensor  Model 

This  section  provides  the  sensor  model  for  proximity 
sensors  that  has  been  incorporated  in  earlier  work  [6,  9]. 
The  received  power  p\ ,  i.e.,  the  power  measurement  of 
the  i-th  sensor  at  time  t,  is  given  by 


Figure  1:  Pd  vs.  the  distance,  where  P/a=0.001,  a  =  2, 
a  =  0.5,  po  =  3000,  ro  =  lm,  and  L  =  100. 


the  distances  to  the  k  targets,  is 


(1) 


where  Pg  k  is  the  power  measured  at  a  reference  distance 
7'o,fc  for  the  fc-th  target,  K  is  the  total  number  of  targets, 
r\  k  is  the  relative  distance  between  the  i-tli  sensor  and 
the  fc-th  target,  a  is  the  attenuation  parameter  that 
depends  on  the  transmission  medium,  and 

v\  ~  N  (fj.v,al)  . 

The  mean  and  variance  of  the  error  v\  is  derived  from 
the  zero  mean  measurement  noise  of  variance  a2  for  the 
case  that  the  measured  power  is  the  result  of  integrating 
the  square  of  the  measurements  over  L  samples.  As 
shown  in  [6,  9],  pv  =  a2  and  a2  =  2 cr4/L.  This  paper 
assumes  the  references  r'o.fc  =  ro  and  po,fc  =  Po  for  all 
K  targets,  and  ro  and  po  are  known. 

The  i-th  sensor  measures  the  received  power  p\,  pro¬ 
cesses  it  locally  and  reports  a  single  binary  digit:  ‘1’ 
for  the  presence  of  one  or  more  targets  or  ’0’  for  the 
absence  of  any  target.  The  decision  follows  the  rule 


d  = 


1  p\  >  A, 
0  p\<  A. 


(2) 


Pd(  r*)  =  Prob(d-  =  l|targets),  (4) 

=  Prob(p-  >  A| targets), 


=  Q 


A  2^k=lPo,k  \rf^J  k-v 


where  r(  = 


rt 

'  i, 0) 


r\  k  .  Fig.  1  illustrates  Pd  as  a 


function  of  the  distance  of  a  single  target  to  a  sensor. 
The  figure  demonstrates  that  Pd  as  given  by  (4)  behaves 
as  a  sigmoid  function  that  steps  down  from  a  value  near 
one  to  a  value  near  Pfa  over  a  transition  region  of  about 
200  meters. 


3  Maximum  Likelihood  Local¬ 
ization 

The  single  target  ML  localization  algorithm  for  binary 
sensors  was  devised  in  [2]  and  further  explored  in  [4]. 
The  likelihood  is  the  probability  of  the  observed  mea¬ 
surements  for  the  N  sensors,  i.e,  d\  for  i  =  1 
conditioned  on  the  location  of  the  targets  Xk  for  k  = 
1 , ,K  and  the  location  of  the  sensor  nodes  5)  for 
i  =  1, . . . ,  N.  Thus,  the  AT-target  likelihood  function  is 
given  by 


Given  a  probability  of  false  alarm  Pfa,  the  threshold 
A  is  computed  as: 

A  =  UvQ~1(Pfa)  +  Hv,  (3) 


N 

l[(Pd([\x1-Si\,...,\XK-Si\])f- 

2=1 


•  (1  —  -Pd  (0®!  —  Si\,.  .  .  ,  \xk  —  <S'j|]))(1_d‘) 


(5) 


where  Q(-)  is  the  Q-function 


Q(x) 


V2n  . 


T dt . 


Furthermore,  the  probability  of  detection  Pd  at  the  i-th 
sensor,  for  a  given  threshold  (or  Pfa ),  as  a  function  of 


One  can  maximize  (5)  to  obtain  the  K  estimated  target 
positions,  which  are  ML  estimates  [7].  For  the  case  that 
K  =  1,  (5)  simplifies  to  the  1-target  likelihood  function 
as  originally  described  in  [2]. 

The  ML  target  location  estimates  are  obtained  by 
using  a  hill  climbing  algorithm  to  maximize  (5).  Be¬ 
cause  the  likelihood  surface  may  have  multiple  peaks, 


it  is  important  to  initialize  the  hill  climbing  algorithm 
with  a  good  initial  guess  for  the  target  locations.  Our 
previous  work  in  [7]  considers  three  practical  initializa¬ 
tion  approaches:  1)  centroid  of  fc-means,  2)  centroid  of 
largest  two  maximal  cliques,  and  3)  peak  in  1-target 
likelihood  surface. 

The  first  two  initialization  methods  employ  cluster¬ 
ing  approaches;  namely,  fc-means  clustering  and  clique 
clustering.  In  general,  both  methods  cluster  the  sen¬ 
sor  detections  in  I\  classes  associated  to  the  known  K 
number  of  targets.  The  first  approach  employs  the  k- 
means  algorithm  by  clustering  based  on  the  geographic 
positions  of  the  nodes  reporting  detections  and  returns 
the  K  centroids  as  the  initial  target  position  estimates. 
On  the  other  hand,  the  maximal  clique  approach  first 
creates  a  graph  of  the  nodes  returning  detections  where 
edges  exist  between  nodes  when  they  are  both  able  to 
detect  a  common  target  at  sufficient  Pd.  Then,  it  finds 
K  maximal  cliques  under  the  assumption  that  each 
clique  is  generated  from  a  different  target.  Finally,  it 
computes  the  centroid  associated  to  the  node  positions 
from  the  K  cliques  as  the  initial  estimates  of  the  target 
position. 

Unlike  the  first  two  methods,  the  1-target  likelihood 
surface  does  not  assume  any  association  between  the 
sensor  detections  and  the  targets.  However,  it  is  com¬ 
putationally  complex.  It  first  assumes  one  target,  i.e. , 
K  =  1  and  calculates  the  one-target  likelihood  of  the 
target  position  over  a  uniform  grid  of  points  by  (5). 
Note  that  the  largest  peaks  on  the  surface  should  corre¬ 
spond  to  the  true  target  locations.  Because  it  is  impos¬ 
sible  to  visualize  the  general  A'-target  likelihood  surface 
when  K  >  1,  the  1-target  likelihood  surface  provides  a 
convenient  means  to  investigate  how  sensor  placement 
can  affect  the  ability  to  resolve  targets.  To  this  end, 
two  targets  are  considered  resolved  when  the  two  peaks 
closest  to  the  ground  truth  target  positions  are  the  two 
largest  peaks  in  the  likelihood  map. 

Fig.  2  shows  the  likelihood  surface  for  two  targets 
separated  by  rs  =  600m.  The  bright  regions  are  ei¬ 
ther  due  to  the  possible  target  locations  or  no  sensor 
coverage,  and  the  dark  blue  disk  regions  are  due  to  sen¬ 
sors  reporting  ’no  targets’.  Figs.  2(a)  and  (c)  show  the 
physical  configuration  of  the  targets  (the  red  pluses), 
the  nodes  that  do  not  detect  targets  (the  blue  dots), 
and  the  nodes  that  do  detect  a  target  (the  larger  green 
faced  circles).  Figs.  2(b)  and  (d)  show  the  likelihood 
surfaces  with  white  pluses  as  cues  for  the  target  loca¬ 
tions.  In  Figs.  2(a)  and  (b),  no  sensors  report  target 
locations  and  small  peaks  occur  in  the  likelihood  surface 
at  locations  beyond  the  detection  range  of  the  sensors. 
In  Figs.  2(c)  and  (d),  two  sensors  do  detect  the  tar¬ 
gets,  but  the  resulting  peaks  in  the  likelihood  surface 
are  modest  because  the  nearby  sensors  that  do  not  de¬ 
tect  the  targets  provide  evidence  against  the  existence 
of  any  targets.  In  either  case,  the  sensor  network  ge¬ 
ometry  is  unable  to  resolve  the  two  targets. 
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Figure  2:  Demonstration  of  the  ability  of  the  1-target 
likelihood  surface  to  resolve  two  targets  separated  by 
600m:  (a)  The  collection  geometry  where  no  sensors 
report  detections,  (b)  the  corresponding  likelihood  sur¬ 
face,  (c)  the  collection  geometry  when  two  closest  sen¬ 
sors  report  detections,  and  (d)  the  corresponding  like¬ 
lihood  surface. 


Fig.  3  shows  the  likelihood  surface  for  two  targets 
that  are  even  closer  than  in  Fig.  2;  namely,  rs  =  300m. 
In  Fig.  3(a)  most  of  the  nodes  within  detection  range 
of  the  targets  do  detect  the  target.  The  corresponding 
likelihood  surface  demonstrates  a  broad  peak  that  cov¬ 
ers  both  targets,  and  the  targets  are  not  resolved.  In 
Fig.  3(c),  a  sensor  that  lies  between  the  two  targets  does 
not  report  a  detection.  In  that  case,  this  “between” 
node  is  able  to  cut  the  peak  in  Fig.  3(b)  into  two  peaks 
in  Fig.  3(d),  and  the  two  targets  are  resolved.  Overall, 
Figs.  2  and  3  demonstrate  that  to  resolve  the  targets, 
most  of  the  sensors  near  the  targets  must  report  detec¬ 
tions  and  a  “between”  sensor  must  report  no  detection. 
This  insight  leads  to  a  proximity  sensor  network  design 
strategy  as  reported  in  the  next  section. 

4  Design  of  Operation  Parame¬ 
ters 

Inspection  of  the  likelihood  surfaces  reveals  the  condi¬ 
tions  necessary  for  the  network  of  proximity  sensors  to 
resolve  two  targets.  Two  targets  are  resolved  when  the 
two  largest  peaks  in  the  likelihood  surface  correspond 
to  the  location  of  the  targets.  Because  it  is  impossi¬ 
ble  to  visualize  the  2-target  likelihood  surface,  we  use 
the  analysis  of  the  1-target  likelihood  surface  to  dictate 
the  conditions  necessary  to  achieve  target  resolution. 
As  revealed  in  Section  3,  two  targets  separated  by  rs 
meters  can  be  resolved  when  1)  at  least  one  sensor  is 
close  enough  to  each  target  to  detect  it,  2)  most  of 
the  sensors  that  can  detect  a  target  actually  report  a 
detection,  and  3)  there  exist  a  sensor  between  the  tar¬ 
gets  that  is  sufficiently  far  from  either  target  to  report 
a  detection.  The  first  two  conditions  can  be  achieved 
by  various  combinations  of  threshold  and  sensor  density 
settings.  By  lowering  the  threshold,  a  smaller  node  den¬ 
sity  can  achieve  the  first  two  conditions.  However,  the 
third  condition  provides  a  lower  bound  for  the  threshold 
so  that  the  “between”  sensor  cannot  detect  target. 

The  section  provides  the  design  criteria  to  obtain  a 
specified  probability  of  resolution  Pr  at  a  specified  tar¬ 
get  separation  distance  rs.  This  criteria  is  for  the  case 
that  the  sensor  nodes  are  uniformly  distributed  over 
its  coverage  area  so  that  the  node  locations  approxi¬ 
mate  a  realization  of  a  Poisson  Point  Process  (PPP) 
[10].  Based  on  the  conditions  of  target  resolution,  the 
probability  of  resolution  Pr  can  be  bounded  by 

Pr  >  (tit2)2(1  -  r3),  (6) 

where  T\  is  the  lower  bound  probability  that  a  sensor 
within  range  r  to  the  target  can  detect  it,  t2  is  the 
probability  that  at  least  one  sensor  is  within  range  r  and 
t 3  is  the  probability  that  the  “between”  sensor  reports 
a  detection.  Because  the  first  two  conditions  must  be 
met  for  both  targets,  t\  and  r2  are  squared  in  (6). 

Actually,  is  an  approximation  for  the  lower 

bound  for  the  first  two  conditions  to  be  met  at  a  single 


Figure  3:  Demonstration  of  the  ability  of  the  1-target 
likelihood  map  to  resolve  two  targets  separated  by 
300m:  (a)  The  collection  geometry  where  13  sensors 
report  detections,  (b)  the  corresponding  likelihood  sur¬ 
face,  (c)  the  collection  geometry  where  12  sensors  report 
detections  while  the  “between”  sensor  reports  ’no  de¬ 
tection’,  and  (d)  the  corresponding  likelihood  surface. 


target.  By  definition, 


n  =  Q 


(7) 


since  by  (4),  Pd  is  larger  for  targets  closer  than  range  r. 
Furthermore,  for  a  given  sensor  density  p  of  the  PPP 
[11,  12] 

t2  =  1  —  exp(-7 Tpr2).  (8) 


effective  detection  range  of  the  sensors.  However,  the 
threshold  cannot  be  made  arbitrarily  small  or  else  the 
final  resolvability  condition  cannot  be  met.  To  meet 
the  final  condition,  the  probability  that  the  “between” 
sensor  reports  a  detection  is 


T3  =  Q 


A- 


(0.5r„)° 


Pv 


(15) 


Note  that  p  =  N/A,  where  N  is  the  number  of  sensor 
nodes  within  the  network  and  A  is  the  area  covered  by 
sensor  network.  Now,  if  we  treat  T\  as  the  probability  of 
detection  for  all  nodes  within  range  r  to  the  target,  then 
the  nodes  that  report  a  detection  and  do  not  report  a 
detection  form  PPPs  whose  densities  are  Tip  and  (1  — 
Ti)p,  respectively  [10].  Thus,  within  the  circular  region 
of  radius  r  around  the  target,  the  probability  that  at 
least  one  node  reports  a  detection  is  given  by 

Pi  =  1  —  exp(— pTiTrr2).  (9) 


Thus,  the  threshold  cannot  be  smaller  than 


\  =  avQ  1(r3)  + 


2porg 

(0.5rs)a 


+  Pvi 


(16) 


and  insertion  of  (16)  into  (14)  leads  to  the  smallest 
achievable  node  density  pc  for  given  values  of  r’s,  i.e., 


Pc  =  -[(Tv(Q  1  (r"3 )  Q  1(n))  + 


(0.5rs)c 


ln(l  —  T2) 


W 


(17) 


Furthermore,  the  probability  that  all  nodes  within  the 
circular  region  report  a  detection  is 

Pa  =  Piexp(-p(l  -  n)7rr2),  (10) 

=  exp(— p{l  —  ti)*k  r2)  —  exp(— pirr2). 

The  probability  that  most  sensors  within  the  circular 
region  detect  the  target  is  bounded  below  by  pa  and 
above  by  pi .  As  shown  in  the  Appendix, 


Pa  <  tit2  <  pi- 


(11) 


As  r 1  increases  to  one,  the  bounds  in  (11)  become  tight. 
Hence,  tit2  is  a  reasonable  proxy  for  the  probability 
that  the  first  two  conditions  of  resolvability  are  met  for 
a  single  target. 

For  a  given  value  of  r2,  the  achievable  node  density 
as  a  function  of  the  detection  range  that  achieves  a  Pd 
of  at  least  n  is  given  by 


P  = 


ln(l  -  r2 ) 

7rr2 


(12) 


where  by  the  inverse  of  (7) 


r  = 


Po  ro 


(A  -  avQ  1(ri) 


(13) 


Thus,  the  relationship  between  achievable  p  and  param¬ 
eters  ri  and  r2  is 


p  =  - 2  T2^  (A  -  avQ  1  (ti )  pv )  ^  .  (14) 

*Po  ro 

To  meet  the  first  two  resolvability  conditions  for  a  given 
r 1  and  r2,  the  node  density  can  be  made  smaller  by 
simply  decreasing  the  threshold  A,  i.e.,  increasing  the 


Judicious  selection  of  the  t's  can  provide  a  small  crit¬ 
ical  density.  For  the  Pr  specification,  (6)  leads  to  two 
degrees  of  freedom  for  the  r’s,  where  by  considering 
arbitrary  n  and  t2,  t3  =  1  —  (T^2p  ■  A  good  design 
should  determine  the  smallest  pc.  For  this  reason,  the 
sensor  network  design  minimizes  (17)  to  obtain  the  crit¬ 
ical  sensor  density: 


*  /  *  *  \ 
Pc  —  Pc  (T1  5  T2  )  ! 


(18) 


where 

{t* ,  r|}  =  arg  min  pc(r 1 ,  r2) .  (19) 

ti,t2 

Inspection  of  (17)  over  a  grid  of  admissible  values  of 
ri  and  r2  indicates  that  there  is  a  unique  minimum. 
Future  work  will  try  to  prove  this  claim.  We  deter¬ 
mine  p*c  using  the  Nelder-Mead  simplex  algorithm  via 
the  fminsearch  function  in  Matlab.  The  search  is  ini¬ 
tialized  by  setting  n  =  t2  =  \fIA-. 

Overall,  to  achieve  a  probability  of  resolution  of  Pr 
for  two  targets  separated  by  a  distance  of  rs  using  a 
proximity  sensor  network  covering  an  area  of  A,  one 
must  use  Nc  =  p*  A  sensors  with  a  threshold  setting  of 


A* 


OvQ 


2porg 

(0.5rs)Q 


+  Pv 


(20) 


Equivalently,  by  (3),  the  P/a  for  the  sensor  node  de¬ 
signed  to  be 


Pfa  =  Q 


Pr 


(n* 


r2*)5 


2 Porg  A 
av(0.5rs)a  J 


(21) 


The  proximity  sensor  network  design  is  accomplished 
via  (17)-(21). 


600 


0 


Figure  4:  Critical  number  of  sensors  versus  specified  rs 
for  various  specified  values  of  Pr  using  (17). 

5  Simulations 

This  section  provides  simulations  when  the  propaga¬ 
tion  loss  factor  is  a  =  2,  the  target  power  is  p0  =  3000 
at  a  normalized  distance  of  rg  =  lm,  the  ambient  noise 
is  cr  =  0.5,  and  the  number  of  integrated  samples  is 
L  =  100.  Furthermore,  the  sensor  field  covers  a  2km  x 
2km  region  so  that  A  =  4km2.  Fig.  4  shows  the  the¬ 
oretical  critical  number  of  sensors,  i.e. ,  the  minimum 
of  (17)  via  fminsearch  multiplied  by  A,  as  a  func¬ 
tion  of  the  specified  separation  distance  rs  for  various 
specified  values  of  Pr.  As  expected,  the  critical  num¬ 
ber  of  sensors  decreases  as  the  separation  distance  in¬ 
creases.  Fig.  5  shows  the  corresponding  critical  Pfa  for 
the  given  design  specifications  as  given  by  (21).  For  a 
given  Pr,  as  the  separation  distance  decreases,  the  Pfa 
must  decrease,  i.e.,  the  threshold  increases,  so  that  the 
’’between”  sensor  does  not  detect  either  target.  For  a 
given  rs,  as  Pr  increases,  more  senors  are  required,  and 
the  Pfa  decreases  slightly  due  to  the  increase  in  sensor 
density.  Note  that  for  rs  <  400,  the  Pfa  is  too  small 
to  be  computed  due  to  the  numerical  precision  of  the 
computer.  However,  the  corresponding  thresholds  A  are 
finite. 

To  prove  that  the  range  of  P/a  and  the  critical  sen¬ 
sor  density  work  to  meet  the  desirable  Pr  at  a  spec¬ 
ified  rs,  we  ran  5000  Monte  Carlo  experiments  com¬ 
posed  of  100  random  sensor  network  configurations  x 
50  Monte  Carlo  runs  per  sensor  network  configuration. 
In  each  configuration,  the  sensors  are  randomly  placed 
in  the  2km  x  2km  field  via  a  uniform  distribution.  The 
number  of  nodes  placed  in  the  network  and  their  cor¬ 
responding  threshold  (or  Pfa)  settings  were  based  on 
the  design  criteria  as  given  by  (17)-(21).  These  values 
are  provided  in  Figs.  4  and  5.  For  each  Monte  Carlo 
realization,  we  compute  the  one-target  likelihood  sur¬ 
face  over  the  whole  region.  In  the  one-target  likelihood 
map,  the  two  targets  are  considered  resolved  when  the 


Figure  5:  Critical  P/a  versus  specified  rs  for  various 
specified  values  of  Pr  using  (21) 


Figure  6:  Percentage  of  targets  resolved  when  the  sen¬ 
sor  density  is  p*  and  the  thresholds  are  set  to  achieve 
a  probability  of  false  alarm  of  p*^a. 

two  peaks  closest  to  the  ground  truth  target  positions 
are  the  two  largest  peaks.  The  results  of  the  actual 
percentage  of  targets  resolved  are  shown  in  Fig.  6.  It  is 
clear  that  the  actual  percentage  of  resolved  targets  ex¬ 
ceed  the  specified  Pr  when  the  targets  are  separated  by 
more  than  600m.  The  designs  are  conservative  as  the 
achievable  probability  of  resolution  should  exceed  the 
specified  Pr  because  of  (6).  When  the  targets  are  sep¬ 
arated  by  400m,  the  actual  Pr  is  lower  than  expected 
when  the  desired  Pr  =  0.8  or  Pr  =  0.9.  This  may 
be  due  to  the  fact  that  the  specification  of  r3  in  Sec¬ 
tion  4  only  provides  the  desired  false  alarm  rate  of  the 
’’between”  sensors  and  does  not  actually  ensure  the  ex¬ 
istence  of  such  a  “between”  node.  At  the  lower  target 
separations  the  area  for  the  ’’between”  sensor  to  exists 
becomes  too  small  to  be  accommodated  by  the  critical 
density  that  achieves  the  first  two  resolvability  condi¬ 
tions. 


Finally,  we  evaluated  the  ML  localization  method 
when  the  sensor  network  is  designed  via  the  resolv¬ 
ability  criteria.  Specifically,  we  ran  the  ML  estimators 
with  various  practical  initialization  methods:  centroid 
of  cliques,  centroid  of  fc-means,  and  peaks  in  the  1- 
target  likelihood  surface.  We  also  ran  the  ML  estimator 
using  the  ground  truth  initialization  as  the  performance 
bound.  Then,  we  determined  the  percentage  of  times 
that  the  practical  initialization  methods  are  able  to  pro¬ 
vide  an  estimate  that  is  comparable  to  ground  truth 
initialization  method.  Fig.  7  shows  the  percentage  of 
times  the  ML  estimates  due  to  the  ground  truth  initial¬ 
ization  agreed  with  the  ML  estimates  due  to  one  of  the 
practical  initialization  methods.  For  most  of  the  time, 
the  correspondence  to  the  ground  truth  initialization 
meets  or  exceeds  the  specified  Pr  when  the  target  sepa¬ 
ration  is  greater  than  500m  for  the  clique  and  1-target 
likelihood  methods.  For  rs  <  500m,  the  ground  truth 
correspondence  drops  similar  to  the  drop  in  achievable 
Pr  as  seen  in  Fig.  6.  When  using  the  fc-means  for  ini¬ 
tialization,  the  ground  truth  correspondence  does  not 
achieve  the  specified  Pr  for  most  cases  when  Pr  >  0.8. 

6  Conclusion 

This  work  documents  our  initial  attempt  to  deter¬ 
mine  the  operating  parameters  of  a  proximity  sensor 
network  to  achieve  a  specified  probability  of  resolution 
at  a  specified  target  separation  rs.  The  probability  of 
resolution  is  derived  by  considering  three  necessary  con¬ 
ditions  to  resolve  two  target.  These  conditions  consist 
of  1)  node  proximity  to  the  targets,  2)  sufficient  number 
of  detections  from  nodes  close  to  the  target,  and  3)  the 
existence  of  the  “between”  node  that  cannot  detect  ei¬ 
ther  target.  These  conditions  lead  to  a  design  strat¬ 
egy  that  determines  the  necessary  sensor  density  and 
threshold  settings  to  achieve  the  desired  Pr  for  a  given 
rs.  Simulations  demonstrate  that  the  design  strategy 
exceeds  the  specified  Pr  when  rs  >  500m.  For  lower 
values  of  rs,  the  realized  Pr  can  be  lower  than  speci¬ 
fied.  This  may  be  due  to  the  fact  that  the  design  does 
not  consider  whether  the  node  density  is  high  enough 
for  the  “between”  sensor  to  exist  with  sufficiently  high 
probability.  In  other  words,  the  design  method  assumes 
that  a  node  at  the  midpoint  of  the  two  targets  cannot 
detect  either  target.  However,  it  does  not  factor  into 
account  the  probability  that  the  node  will  actually  lie 
near  this  midpoint. 

Future  work  will  revise  the  design  criteria  to  still 
achieve  the  desired  Pr  when  the  target  separation  is 
small.  This  may  be  accomplished  by  factoring  in  the 
probability  that  a  “between”  sensor  node  exist.  Fur¬ 
thermore,  we  plan  to  investigate  methods  to  predict  the 
root  mean  squared  position  error  performance  of  the 
ML  localization  estimate.  Finally,  the  simplifying  as¬ 
sumptions  can  be  relaxed  to  obtain  performance  bounds 
under  ever  increasing  realistic  conditions. 


Separation  Distance  (m) 


(c) 

Figure  7:  Percentage  of  times  the  practical  initializa¬ 
tion  leads  to  an  ML  estimate  that  corresponds  to  the 
ML  estimate  when  initializing  by  the  ground  truth:  (a) 
Clique  initialization,  (b)  fc-means  initialization,  and  (c) 
1-target  likelihood  initialization. 


A  Proof  for  the  Inequalities  in 

(H) 

Let  a  =  exp(— pirr2).  Note  that  a  is  bounded  by 
0  <  a  <  1.  Then  by  considering  (8),  the  bounds  in  (11) 
can  be  expressed  as 

a1-”1  —  a  <  n(l  —  a)  <  1  —  aTl.  (22) 


Since  n  corresponds  to  the  probability  of  detection,  it’s 
bounded  by  0  <  n  <  1.  Let’s  define  a  function  of  n  as 

fin)  =  a1_Tl  -  a  -  n(l  -  a).  (23) 


Clearly  /( 0)  =  /(l)  =  0.  The  second  derivative  of  /(•) 
with  respect  to  T\  is 

=  (ln(a))2a1_ri  >  0.  (24) 

Thus,  the  function  is  convex,  which  means  that 

fir, i)  <  (1  -Tl)m  +  nf(l),  (25) 

<  0, 


and 

a 1~ri  —  a  <  ri(l  —  a).  (26) 

This  proves  the  first  inequality  in  (22).  Similarly,  we 
can  define  a  function  of  t\  as 


fin)  =  ri(l  -  a)  -  1  +  aTl. 


(27) 


Clearly,  /( 0)  =  /( 1)  =  0  Furthermore  the  second 
derivative  of  /(•)  with  respect  to  n  is 


H 

d2Ti 


=  (ln(a))2a[  >  0. 


Again  the  function  is  convex,  and 

fin)  <  (1  -  n)/(0)  +  ti/(1), 

<  0. 

Thus, 

ri(l  —  a)  <  1  —  aTl, 

which  proves  the  second  inequality  in  (22). 


(28) 


(29) 


(30) 
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